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Colossal enhancement of upper critical fields in noncentrosymmetric heavy fermion 
superconductors near quantum criticality: CeRhSi3 and CeIrSi3 

Y. TadaB N. Kawakami, and S. Fujimoto 

Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

Strong coupling effects on the upper critical fields H C 2 along the c-axis in the noncentrosymmetric 
heavy fermion superconductors near quantum criticality CeRhSi3 and CeIrSi3 are examined. For 
sufficiently large spin-orbit interactions due to the lack of inversion symmetry, H C 2 is mainly de- 
termined by the orbital depairing effects. From microscopic calculations taking into account the 
strong spin fluctuations, we show that H C 2 is extremely enhanced as the system approaches the 
quantum critical point, resulting in H C 2 ~ 30T even in the case with the low transition temperature 
T C (H = 0) ~ IK, which well explains the huge H c % observed in the recent experiments. Our results 
reveal intrinsic and universal properties common to strong coupling superconductivity caused by 
spin fluctuations near quantum criticality. 

PACS numbers: 



Quantum criticality is a phenomenon emerging in the 
vicinity of a second-order phase transition at sufficiently 
low temperatures governed not by thermal fluctuations 
but by quantum fluctuations. There has been accu- 
mulated evidence that quantum criticality plays an es- 
sential role in unconventional superconductivity realized 
in strongly correlated electron systems such as heavy 
fermion compounds and high-T c cuprates 1 . In partic- 
ular, it has been proposed by several authors that mag- 
netic critical fluctuations mediate Cooper pairing in some 
classes of unconventional superconductors2^' 4 *£&£. In 
this scenario, it is expected that the pairing interac- 
tion is explosively enhanced just at the quantum critical 
point (QCP), though such enhancement has not been ob- 
served so far in real systems, partly because the QCP at 
zero temperature associated with the magnetic order is 
veiled by the superconducting phase. In this respect, the 
recent experimental studies on the upper critical fields 
of noncentrosymmetric heavy fermion superconductors 
CeRhSi^^ and CelrSi j 10 ' 11 under applied pressure are 
very curiou o 12 ! 13 . In these systems, the distance from 
the QCP is controlled by the applied pressure. The re- 
markable features of the experimental results are as fol- 
lows: (i) As the pressure approaches a critical value, these 
systems exhibit extremely high upper critical fields H C 2 
which exceed the orbital limit as well as the Pauli limit 
estimated by the conventional BCS theory. The observed 
H C 2 ~ 30T is, surprisingly, almost comparable to those 
of high-T c systems, though the transition temperature at 
zero field is merely ~ IK. (ii) The increase in the upper 
critical field is enormously accelerated as the temperature 
T is decreased, making a sharp contrast to any other 
superconductors in which the increase of H C 2 becomes 
slower as T is decreased, (iii) H C 2 increases very rapidly 
as the pressure approaches the critical value, though, by 
contrast, the pressure dependence of T c is moderate. 

The absence of the Pauli depairing effect in these 
systems is consistent with the previous theoretical 
predictions 14 . However, the origins of the strong suppres- 



sion of the orbital depairing effect and the other extraor- 
dinary behaviors (ii) and (iii) have not been elucidated 
so far. It is particularly important to clarify how these 
unusual properties are related to quantum criticality and 
also what kind of role is played by the noncentrosym- 
metry in the compounds CeRh(Ir)Si3. The answers to 
the questions should certainly provide us with deep un- 
derstanding of superconductivity for strongly correlated 
electrons near quantum criticality. 

In this letter, we demonstrate that the above striking 
features are naturally understood in terms of the strongly 
enhanced pairing interaction due to spin fluctuations in 
the vicinity of the magnetic QCP. Our results reveal that 
the above features experimentally observed for CeRhSi3 
and CeIrSi3 are not specific to these noncentrosymmet- 
ric superconductors, but rather intrinsic and universal 
in strong coupling superconductivity caused by critical 
fluctuations, as long as the Pauli depairing effect is sup- 
pressed by some other mechanisms such as the lack of 
inversion center of the crystal structure as in the case of 
CeRh(Ir)Si3. Therefore, the unusual properties found for 
the upper critical fields uniquely characterize the inter- 
play between quantum criticality and spin-fluctuation- 
mcdiatcd Cooper pairings. 

We focus on the low energy quasiparticles in the clean 
limit mainly formed by /-electrons and approximate 
the system by the following effective single-band tight- 
binding model, 

S = So + Ssf, 

S = 4 (-iu n + e(fc)) Ck + ^ c' k aCo(k) ■ crc k , 

k k 
kk' q 

where Cfe = (cfe-[-,Cfej_) is the annihilation operator of 
the Kramers doublet. Here we have introduced the 
notation k = (iuj n ,k). The lack of inversion sym- 
metry leads to the spin-orbit interaction which is rep- 
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resented in the second term with the Zeeman split- 
ting. The third term is the electron-electron inter- 
actions through spin fluctuations. For CeRhSi3 and 
CeIrSi3 which have body centered tetragonal lattice 
structures, the dispersion relation and the spin-orbit in- 
teraction are given by e(k) = —2ti(cosk x + cosk y ) + 
4i2 cos k x cos k y — cos(k x /2) cos(k y /2) cos k z — [i and 
oiCq = (asinfcy, — asinfc^, ^^bH). We fix the parame- 
ters as (tx,t2,t 3 ,n,a) = (1.0,0.475,0.3, 1.05,0.5) by tak- 
ing t\ as the energy unit. The Fermi surface deter- 
mined by these parameters is in qualitative agreement 
with the band calculation and can reproduce the peak 
structures of the momentum-dependent susceptibility ob- 
served by the neutron scattering experiment s 16 ' 17 ! 18 . To 
describe the strong spin fluctuations near the QCP, we 
phenomenologically introduce the effective interaction 
between quasiparticle s 3 ' 19 ' 20 ' 21 , 

x{Wn ' q) = ?i + e 2 (q-QJ 2 + KI/(r r 2 ) 

and £(T) = ^£+8 ' wnere Xo and Tq are respectively 
the susceptibility and the energy scale of spin fluctua- 
tions without strong correlations. These quantities are 
renormalized through the coherence length £(T) as the 
system approaches the QCP. The critical exponent of £ 
is the mean field value 1/2 and 9 is considered to de- 
crease monotonically as the applied pressure approaches 
the critical value for the AF orde r 19 ' 20 . The tempera- 
ture dependence of £ is also consistent with the recent 
NMR experiment for CelrSijr 2 ^,. The propagating vec- 
tors are Q l = (±0.43tt, 0, 0.5tt), Q 2 = (0, ±0.43tt, 0.5tt) 
according to the neutron scattering experiments^ 7 -. In 
this study, we fix the parameters in x(iv ni q) as Tq = 3.6 
and £ = 0.4647. The former is of the same order as the 
Fermi energy and the latter is determined so that the 
maximum of £ would be £ m ax ~ Sk^ 1 , where 5k is the 
interstice of the fc-mesh of the Brillouin zone used in our 
numerical calculations. 

To describe the quasiparticles in the strong cou- 
pling regime, we introduce the normal selfenergy up to 
the first order in g 2 Xo which is given by S SlS2 (fc) = 
Ji Efe' 9 2 x(k ~ k')G° SlS2 (k')6 SlS2 , neglecting off diagonal 
elements which are not important in the present study. 
G° is the Green's function for noninteracting quasiparti- 
cles with spin indices s\, S2- ReE only gives the deforma- 
tion of the Fermi surface while ImE gives two important 
effects to the quasiparticles around the Fermi level; the 
mass enhancement by the factor z^ 1 = 1 — (ImE(7rT) — 
ImE(— ttT))/2ttT and the damping factor 7 = —ImE. It 
is reasonable to consider that e(k) already includes the 
shift due to ReE and to replace e(fe) + ReE(fc) — > e(k). 
We also neglect the change in Cq by the normal selfen- 
ergy, which has little influence in the discussion of the 
upper critical fields. We note that the coupling constant 
g should be regarded as an effective one renormalized by 
the vertex correction s 23 ' 24 . 

To study the transition temperatures under applied 



fields, we use the linear lized Eliashberg equations where 
the Green's functions of the quasiparticles depend on the 
vector potential A. We employ the familiar quasiclassical 
approximation for the Green's function G(iuj n , x, y; A) = 
e ie*(x,V) G (iu} n ,x-y;A = 0), <p(x,y) = J* A(s)ds 
which is legitimate for kpln S> 1 where fcp is the Fermi 
wave number and Ih the magnetic length. We take 
only the N = Landau level because, as shown below, 
H C 2 is mainly determined by the orbital depairing ef- 
fect. The gap function is approximated by A aa i(k; R) = 
A aQ ' (k)(po(R x , Ry), where <fio is the lowest Landau level 
wave function. The resulting Eliashberg equations for 
A(fc) with applied field H = (0, 0, H) are, 

A aa ,(k) = — ^2 Vaa',p/3'(k,k') 

fc',/3/3'77' 

/ 1 + rC"{k) • er\ 

}^ I 2 J * sgn ^ ^ 

(^) 1/2/ (OT) G ^' ( - fc ' )A -' (fc ' )(1) 

where \k) = C%{k)/ y/^Li ££(*0 2 > = w„ - 

ImE(fc), a T (k) = ^\H{v TX {k) 2 + v Ty {k) 2 ), b T (k) = 
\u)(k)\ +isgn(ui)e T (k), with the dispersion e T (k) — e(k) + 

tol \JY^h=i £-oife) 2 anc ^ tne velocity v T (k) = Ve T (fe). 

f(z) is defined as f(z) — erfc(z). Here, we have 

neglected the II operators in 7 = — ImE because V7 • II 
represents the contribution to H C 2 from the anisotropy 
of the lifetime of the quasiparticles on the Fermi surface 
and is of secondary importance compared to 7 itself. The 
velocity can be renormalized by the factor z(k), leading 
to enhancement of H C 2 . The effective pairing interaction 
V is evaluated at the lowest order in g 2 Xo, 

V ss Mk,k') = -\g 2 x{k-k') + ^g 2 x(k + k'), 

V sS Mk,k') = lg 2 x(k-k') + ^g 2 x(k + k'), 

Vss,ss(k,k') — —Vss^sikjk') and the other components 
are zero. The applied fields will not affect V, because 
Tq = 3.6 is large enough compared to the energy scale of 
the applied fields. 

By solving eq. (JTJ) at H = 0, we find that the sym- 
metric superconducting state is most stable among the 
five irreducible representations of point group , which 
is consistent with the previous studyi£. In the follow- 
ing, we focus on this state. The Ai symmetric order 
parameter is in the form of A(icj„, fc) = [A. s (iu> n )do(k) + 
A t (iu> n )d(k) ■ <r]i<j2 with do = cos(2fc z ), d = cos(2fc z )£ - 
In our model, however, the amplitude of the triplet part 
A t is so small (< 0.01A S ) that we will neglect it in the 
following. In FiglU we show the transition temperatures 
T c for the Ai symmetric superconducting state at H = 
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FIG. 1: Transition temperatures for the extended s-wave state 
as functions of g 2 Xo f° r several 9. 



FIG. 3: H% b versus T calculated with T-dependent £ (solid) 
and T-independent £ (dotted) for 9 = 0.002, 0.03. 
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FIG. 2: Temperature dependence of the upper critical fields 
for the extended s-wave state. The dotted line with open 
circles represents the Pauli limiting field for 9 = 0.002 and 
the others the orbital limiting fields for several 9. The unit of 
9 is tx = 113K. 



as functions of g 2 Xo f° r several 9. The transition tem- 
peratures saturate for large g 2 Xo because of the strong 
coupling effects. Their dependence on 9 is quite weak. 

We now turn to the case with H ^ 0. In Fig[5J we 
show upper critical fields as functions of temperature for 
several 9 at g 2 Xo == 13- For this value of g 2 Xo, the renor- 
malization factor z -1 |t=o.oi is ~ 2, which indicates that 
the system is in the strong coupling regime. In Fig. 2, 
we use the following values of the parameters; (lattice 
constant) =4Ap^iii and t\ — 113K. The latter is obtained 
by identifying the calculated maximum transition tem- 
perature T c = 0.0115 as 1.3K which is the averaged value 
of maximum T c for CeRhSi^ and CelrSi^. The dotted 
line with open circles in Fig. 2 corresponds to the Pauli 
limiting field for 9jt\ = 0.002 and the other lines to the 
orbital limiting fields for several 9. In the present param- 
eters, transition temperatures are typically T c ~ O.Olii 
and thus the spin-orbit coupling constant is very large 
compared to T c (a = 0M C ~ 50T C ), leading to ff c p 2 auli 
largely exceeding H° 2 h . Note that, with sufficiently large 
a, -ff^ aull (0) calculated from the Eliashberg equations 
coincides with the one obtained from thermodynamic 
considerations 1 -. From the result, we can conclude that, 



for sufficiently large a>T c , the upper critical fields are 
determined by the orbital limiting fields in accordance 
with the previous studie a 14 ' 15 . A remarkable point in 
Figj2]is that the orbital limiting fields are enhanced as 9 
is decreased, and reach 30T for sufficiently small 9 even 
when T C (H = 0) ~ IK, which well explains the exper- 
imentally observed huge H c2 for CeRhSi3 and CeIrSi3. 
Also, the H c2 curve has a upward curvature, which is 
consistent with the experimental observations. These re- 
sults verify that the superconductivity in these systems 
is mediated by the strong antiferromagnetic spin fluctu- 
ations near the QCP. 

The origin of the colossal enhancement in H° 2 h is un- 
derstood as follows. As the temperature approaches ab- 
solute zero, the antiferromagnetic correlation length £(T) 
increases, leading to the enhancement of the effective 
pairing interaction V oc £(T) 2 , while the pair breaking 
effect due to the quasiparticle damping ImE(T) is sup- 
pressed. As a result, the coherence length of the super- 
conducting state is extremely reduced at low tempera- 
tures. This behavior is significant especially for small 9, 
giving rise to the huge values of H° 2 h , as seen in Fig|2] 
To demonstrate the importance of the rapid increase in 
£(T), we also calculate H c2 by assuming that £(T) is 
independent of temperature and has a constant value, 
£ = £(T C ) (Fig. E}. We see that the H°^° curves cal- 
culated with temperature independent £ exhibit neither 
significant upturn nor enhanced values of H C 2(T) at low 
temperatures. The comparison of these results with those 
computed with the correct £(T), as shown in Fig. [3j con- 
cludes that the rapid increase in £(T) is essential to ob- 
tain the huge values of i/°| b as well as the characteristic 
upturn profile at low temperatures. 

Another remarkable feature of the experimental ob- 
servations for CeRhSi3 and CeIrSi3 is the strong pres- 
sure dependence of H c 2 12 ^ 13 ; i.e. H c2 increases explo- 
sively as the pressure approaches the critical value, while 
the pressure dependence of T c is moderate. To discuss 
this point in detail, we focus on T C (H = 0) and H c2 at 
low temperature as functions of 9. Figure [4] shows the 
9 dependence of the transition temperatures at H = 
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FIG. 4: 9 dependence of the transition temperature t c (9) and 
the upper critical field h C 2(0). The definitions of t c (8) and 
h C 2{9) are given in the text. 

and the upper critical fields at T = T m in = 0.001 nor- 
malized by the values calculated with 6 = (9 max = 0.03; 
t c (6) = T C (H = 0,6)/T c (H = 0,8 = max ) and h c2 {8) = 
H c2 (T = T min ,9)/H c2 {T = T min ,8 = 8 max ). We see that 
t c (6) is slightly changed as 8 is decreased. By contrast, 
h C 2(8) is strongly enhanced, which is consistent with the 
above-mentioned experimental result o 12 i 13 . These behav- 
iors are understood as a result of the strongly enhanced 
pairing interaction at low temperatures in the vicinity of 
the QCP. Although the increase of T C (H = 0) near the 
QCP (i.e. 8 = 0) is considerably suppressed by the pair 
breaking effects due to the inelastic scatterings with spin 
fluctuations, H c2 at low temperatures is not seriously af- 
fected by them, resulting in the strong enhancement of 
H c2 for 6»-0. 

We have theoretically clarified the origins of sev- 
eral puzzling phenomena posed experimentally for 
CeRh(Ir)Si3. Here we wish to stress that the extraor- 
dinary properties of H c2 discussed here are not specific 
to these compounds, but are inherent in the orbital lim- 
ited superconductors caused by spin fluctuations in the 
vicinity of the QCP. However, such anomalous behav- 
iors of H c2 have not been observed so far in any other 



systems than CeRh(Ir)Si3, although there are many su- 
perconducting systems which are believed to be located 
in the vicinity of a magnetic QCP. The key to resolve 
this apparent contradiction is the lack of inversion sym- 
metry. In the presence of the sufficiently large anisotropic 
spin-orbit interaction, as discussed above, the dominant 
mechanism which determines the properties of H c2 is the 
orbital depairing. In addition, generally, the orbital lim- 
iting fields are much more influenced by the electron cor- 
relations than the Pauli ones. In ordinary centrosym- 
metric superconductors, the enhanced strong coupling ef- 
fects on H c2 near a QCP, which could exist in principle, 
are masked by the Pauli depairing effects. The noncen- 
trosymmetric heavy fermion superconductors CeRhSi3 
and CeIrSi3 with large anisotropic spin-orbit interactions 
are the first systems in which H c2 directly exhibits the 
enhanced strong coupling effects due to spin fluctuations 
near the QCP. We thus naturally expect that the anoma- 
lous properties of H c2 could be more generally observed 
in other strong-coupling superconductors caused by crit- 
ical spin fluctuations if the Pauli depairing effect could 
be suppressed by some other mechanisms. 

In summary, we have studied the strong coupling ef- 
fects on the upper critical fields in CeRhSi3 and CeIrSi3. 
We have found that the upper critical fields are extremely 
enhanced by critical spin fluctuations near the QCP, 
which successfully explains the recent experimental ob- 
servations. We have elucidated that the huge values of 
H c2 and its strong dependence on the applied pressure 
is a universal property inherent in the orbital limited su- 
perconductors near the QCP. 
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